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Multidirectional  Functionally  Graded  Plates  Exposed  to  Transient-Type 

Pressure  Pulses 
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ABSTRACT 

Multidirectional  functionally  graded  simply  supported  plates  where  the  grading  occurs 
symmetrically  in  all  three  coordinate  directions  exposed  to  time-dependent  pressure  pulses  is 
considered.  Within  this  context,  the  classical  thin  plate  theory  is  adopted  (CPT).  The  grading 
occurs  in  all  three  coordinate  directions  and  is  based  on  3  independent  polynomial  and  power 
law  distributions.  The  equations  of  motion  are  derived  through  the  use  of  Hamilton’s  Principle. 
By  leveraging  both  approximate  analytical  and  numerical  techniques  such  as  the  Guassian- 
Quadrature  method.  The  Galerkin-Method,  and  the  4*  Order  Runge-Kutta  method,  the  dynamic 
response  for  various  kinds  of  pressure  pulses  is  presented.  In  a  detailed  fashion,  the  influence  of 
the  volume  fraction  indexes,  various  geometrical  and  material  parameters,  and  the  damping  on 
the  dynamic  response  is  presented  and  analyzed.  Finally  validations  are  made  with  specialized 
cases  found  within  the  literature. 
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1.  Introduction 

Plate  and  shell-type  structural  elements  find  themselves  in  all  kinds  of  structures  such  as 
buildings,  automobiles,  and  aircraft/spacecraft.  Their  structural  design  should  be  based  on 
considerations  such  as  temperature,  material  composition,  geometrical  attributes,  and  the  types  of 
loadings  that  they  will  be  exposed  to,  whether  static  or  dynamic.  There  is  an  exhaustive  literature 
base  devoted  to  the  structural  response  of  isotropic,  laminated  composites,  unidirectional  and  bi¬ 
directional  functionally  graded  plates  and  shells  under  various  kinds  of  loadings  based  on  the 
classical  plate  theories  to  the  most  advanced  plate  and  shell  theories.  Although  there  remain  gaps 
that  need  to  be  filled  within  the  present  literature  base,  a  new  frontier  that  yet  remains  to  be 
explored  and  developed  is  the  idea  of  multi-directional  functionally  graded  plates  where  the 
material  properties  vary  in  all  three  coordinate  directions  as  opposed  to  unidirectional  or 
bidirectional  orientations.  This  brings  to  the  forefront  the  idea  of  spatial  tailoring;  whereby  the 
distribution  of  the  constituent  materials  in  all  three  coordinate  directions  can  be  manipulated  to 
achieve  an  optimal  material  distribution  to  achieve  an  optimum  structural  design  to  maintain  its 
structural  integrity  within  its  operating  environment.  This  concept  is  analogous  to  structural 
tailoring  for  laminated  unidirectional  composites  whereby  the  ply-angle  is  manipulated  to 
achieve  an  optimum  configuration. 

Based  on  a  comprehensive  literature  review,  a  theoretical  base  which  to  build  upon  is  found  to 
be  nonexistent.  A  brief  synopsis  on  how  to  approach  and  theoretically  model  multidirectional 
functionally  graded  plates  incorporating  thermal  effects,  with  the  idea  of  spatial  tailoring,  has 
been  introduced  by  [2]  and  [3] .  With  this  in  mind,  a  basic  theoretical  foundation  of  the 
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multidirectional  functionally  graded  plates  based  on  a  linear  classical  thin  plate  theory  under 
various  types  of  dynamic  loadings  for  the  case  of  simply  supported  plates  is  developed.  Within 
this  theoretical  foundation  are  the  inclusion  of  damping  effects  and  the  transverse  inertia. 


2.  Basic  Assumptions  and  Preliminaries 

Shown,  in  Fig.  1,  is  a  pictorial  representation  of  a  multidirectional  functionally  graded  plate 
referred  to  an  Orthogonal  Cartesian  Coordinate  System  (x,  >>,  z) ,  where  z  is  measured  positive 
in  the  upwards  direction  from  the  mid-surface  of  the  plate.  While,  h  is  the  uniform  thickness  of 
the  plate.  Let  any  two  constituent  materials  comprise  a  functionally  graded  plate.  Then  a  generic 
property  P  (x,  >>,  z)  can  be  expressed  in  terms  of  the  volume  fraction  as 

P(x,y,z)  =  [Pi-P2Wi(x,y,z)  +  P2-  (1) 

Where, 

Vi(x,y,z)  =  V^^ix)V^iy)V^^(z).  (2) 


The  chosen  functional  relationships  for  the  volume  fractions,  (x),  (y),  and  (z)  are 


given  in  a  polynomial  and  power  form  as 
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(3c) 


(4) 


2  ’  ^4  M  are  referred  to  as  the  volume  fraction  indexes  providing  a  measure  of  the 
variation  of  the  material  profile  through  the  structure  in  all  three  coordinate  directions.  It  should 
be  noted  that  this  chosen  grading  of  the  constituent  materials  throughout  the  plate  expresses  a 
symmetric  distribution  in  all  three  coordinate  directions.  Also,  depending  on  the  3D  grading 
desired,  other  possible  functional  relationships  for  the  volume  fraction  are  possible  such  as  in  the 
anti-symmetric  case. 

A  detailed  graphical  depiction  of  each  of  the  three  directional  volume  fractions  given  as 
function  of  the  volume  fraction  index  and  their  respective  coordinate  direction  is  provided  in  Fig. 
(2)-(4)  below.  Looking  at  Fig.  (2),  which  displays  the  volume  fraction  as  a  function  of  the 
distance  along  the  plate  in  the  x-direction  for  various  magnitudes  of  the  volume  fraction  indexes 
,  it  is  apparent  that  as  Aj  ^  0  ,  the  volume  fraction  measure,  V^^(x)  ^  1  .  As  Aj  ^  oo  or 
becomes  much  larger  than  zero,  (x)  ^  0  ;  While  other  volume  fraction  indexes  produce 

volume  fractions  falling  between  these  two  extremes.  The  same  exact  behavior  is  seen  in  Fig.  (3) 
where  the  volume  fraction  is  presented  as  a  function  of  the  distance  along  the  plate  in  the  y- 
coordinate  direction  for  various  volume  fraction  indexes,  A  2  .  In  Fig.  (4),  The  volume  fraction  as 
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a  function  of  the  through-the-thickness  coordinate  in  the  z-coordinate  direction  for  various 
volume  fraction  indexes  is  depicted.  When  M=0,  the  volume  fraction,  (z)  =  1 ;  While 

(z)  =  0  as  M  ^  00  or  becomes  very  large.  All  other  volume  fraction  index  values  for  M,  give 

volume  fractions  in  between  the  two  extremes. 

With  this  information  in  hand,  Table  1.  below,  gives  a  summary  for  the  possible  types  of 
functional  grading,  depending  on  the  measure  of  the  volume  fraction  index. 

The  chosen  constituent  materials  for  this  paper  are  ceramic  and  metal.  This  leads  to  the 
expression  of  the  material  properties  given  by 

{E{x,  y,z),  p{x,y,zy\  =  ’  Pan  ]Ti(x,  V,  z)  +  [£„,  ]  (5) 

Where, 

^  cm  ~  ^  c  ^  m ’’  P  cm  ~  P  c  Pm' 


The  variation  of  Poisson’s  ratio,  v{x,  y,  z) ,  is  approximated  as  being  constant  throughout  the 
material  grading  of  the  structure. 

3.  The  Strain-Displacement  Relationships 

For  the  classical  plate  theory,  the  strain-displacement  relationships  are  provided  as  [6]  (7) 
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yxz  =  yyz  =  (8d) 

Where,  (uq,vq,wq)  are  2-D  displacement  quantities  of  the  mid-surface  of  the  plate. 


4.  Constitutive  Equations 

The  constitutive  equations  for  a  point-wise  isotropic  material  for  plane  stress  are  given  by 


^xx 
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Where  the  material  stiffnesses,  Q^-  (x,  y,  z) ,  (i,  J=l,2,6)  are  given  by 


^  ^  £(.r,v,z)  ^  vE(x,y,z)  ^  E(x,y,z} 

Cll  -  C22  -  ^  ’  C12  -  ^  ’  C66  - 


1- V' 


1- V' 


2(1 +  v) 


(10) 


5.  Equations  of  Motion 

Adopting  an  energy  approach,  the  equations  of  motion  are  derived  through  the  use  of  Hamilton’s 
Principle  which  is  expressed  as 


(5/ 


(T  -U  +V)dt 


=  0 


(11) 


Where  are  two  arbitrary  instants  in  time.  U  denotes  the  strain  energy,  V  denotes  the  work, 
and  T  denotes  the  kinetic  energy  of  the  structure,  while  <5  is  the  variational  operator.  Without 
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explicitly  providing  the  individual  energy  terms  they  are  assumed  into  the  energy  functional 


where  the  energy  functional  is  expressed  as 


h 

t.  — 

+<^yy^^yy  +  xyT  xy^^zd  D.  +  f  p  {x ,  y)w  q5w  ^dD. 

-- 

"  L  2 

r  1 

+  J  [Pf{x,  y ,t)  -  p{x,  y)w q)]5w QdD.  \dt  =  0 


(12) 


Q  denotes  the  mid-surface  area  of  the  plate,  P,  is  the  distributed  force  at  the  top  surface,  and 
p{x,  y)  is  the  damping  coefficient  per  unit  area  of  the  plate. 

Considering  Eqs.  (8a-d),  (9),  and  (12),  and  carrying  out  the  integration  throughout  the 
thickness,  integrating  by  parts  where  ever  feasible,  and  taking  into  consideration  the  arbitrary 
and  independent  character  of  variations  results  in  three  equations  of  motion  along  with  the 
associated  boundary  conditions  in  terms  of  displacements.  The  equations  of  motion  are 
determined  as 
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4  dD-^i{x,y)  diy^  Wq)  (.r,  y)  5(V  )  2  2 

D^^{x,y)y  Wq  +  2 - - - - y-+  2 - - - - —+VD^^{x,y){V  wQ) 


dx  dx  dy  dy 

[  D]^]^{x,y)  d^WQ  ^  (x,  j)  {x,  y)  d^WQ  ^ 


-d-v) 

1^  dx^  dy‘ 

+  /j{x,y)wQ  =  Pf{x,y,t) 


dxdy  dxdy 


dy- 


dx" 


+  P(x,  y)wQ 


(15) 


As  a  result  of  the  decoupling  of  the  first  two  equations  of  motion  which  are  in  terms  of  the  in¬ 
plane  displacements,  the  third  equation  of  motion  in  terms  of  the  transverse  displacements 
becomes  the  governing  equation  of  motion.  The  reduced  governing  system  of  equations  along 
with  the  associated  boundary  conditions  becomes 


^  ,  _4  ^dDy^ix,y)  div\o)  ^dD^^ix,y)  diV^o) 

owq  Dii(x,y)V  wq+2 - h  2 - h 


dx 


dx 
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d^Djl  (x,  y)  d^WQ  ^ 
dy^  dx^ 


+  P(x,  y)wQ  +  iu{x,  y)wQ  =  P^{x,  y,t)  . 


For  simply  supported  boundary  conditions,  the  boundary  constraints  are  expressed  as 
Along  the  edges  x  =  (0,  Lj ) 


2  2 

d  Wn  d  Wn 
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dx  dy 


(17) 


Along  the  edges  x  =  (0,  L^) 
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Wq  =  0,  -  £>22  (x,  y) 


O  Wr 


dy- 


-  D^2  3^)  ■ 


O  W( 


=  0. 


(18) 


The  global  stiffness  quantities  appearing  within  the  equations  of  motion  are  defined  as 

h 

[A-j  (x,  y),  B-j  (x,  y),  D-  {x,  y)]  =  j  2^  gy  {x,  y ,  z){\,  Z,  )dz  ,  {i,  j  =  1,  2,  6)  (19) 

2 

For  the  case  of  symmetric  grading  in  all  three  coordinate  directions,  ]  =  0  .  With  the  use  of 

Eqs.  (5)  and  (10),  and  carrying  out  the  indicated  integration  within  Eq.  (19),  the  global  stiffness 
quantities  can  be  expressed  as 


(Ay  ,  By  ,  Z)y  )  =  - ,  ^2  ^  ^3  )  j  = 


(^66^B^,D^)=  /  £3) 
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Where, 
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12 


(20a) 

(20b) 


(21a) 

(21b) 

(21c) 


6.  Solution  Methodology 

As  part  of  the  solution  process  of  Eq.  (16),  the  Galerkin  Method  is  leveraged  which  requires 
that  both  the  essential  (kinematic)  and  natural  boundary  conditions  be  fulfilled.  To  achieve  this 
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requirement,  the  transverse  displacement,  wq  (x,  y)  and  the  transverse  pressure,  Pj  (x,  y,z)  are 
both  expressed  in  terms  of  a  Navier-type  expression  given  as 

w{x,y,t)  =  W  iu„y)  (22) 

Where  2^  =  niK  jLy,  =  n7rlL2  and  (m  ,n)  are  the  number  of  sine  half-waves  in  the 
corresponding  directions;  While, 

P^{x,  y,t)  =  (f)sin(  A^x)sm{  ju^y)  (23) 


Where, 


4  c^2 

=  -  P^{x,  y,  t)  X^x)fim{  /u^y)dxdy 

L1L2  •'0  •'0 


(24) 


with  Pj  {x,  y,  t)  =  Pj  (r) ,  and  integrating  gives 


=  (25) 

Where,  P^  (0  is  the  desired  transverse  pressure.  With  these  expressions  in  hand  and  applying  the 
Galerkin  Method  to  Eqs  (16),  utilizing  Eqs.  (22),  (23),  and  (25),  gives 

16P,(0 

^ mn^ mn  (G  +  P-  mn^ mn  (^)  ^ mn^ mn  (^)  ~  2  ’  (26) 

mn  n 

Where  (f)  is  the  amplitude  of  deflection,  is  the  plate  inertia,  is  the  overall 
damping  coefficient  of  the  plate,  and  K  is  the  overall  stiffness  of  the  plate  all  of  which  are 
given  as 
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Within  the  above  expressions, 


=  cos,{  m  7VX !  L^)  ,  =  sm(  m;rx/Lj) 


(30a, b) 


c„y  =  cos 


in7rylL2),  s  =  sini  nny / L2) 


(31a,b) 


Normalizing  Eq.  (26)  results  in. 


^ mn  (^)  2A  (W  it)  +  (^)  ^mn  (^) 


where  is  the  natural  frequency  of  the  plate,  =  C^„  /2/ 


mn  ®  mn  ^hc 


overall  normalized  damping  coefficient,  and  (f)  =  16  (f) / mn  I is  the  normalized 

transverse  pressure. 

7.  Transverse  Pressure  Loadings 
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In  a  dynamic  world,  structures  can  experience  many  different  types  of  external  excitations  from 
a  variety  of  causes  or  sources.  Whether  constant  or  time- varying  they  are  the  driving  force  for 
their  dynamic  motion.  In  the  supersonic  and  hypersonic  world,  aircraft  can  experience  various 
types  of  pressure  waves  which  can  appear  as  violent/hostile  disturbances  which  can  materialize 
from  a  variety  of  environmental  conditions.  Sonic  boom  is  one  such  example  which  occurs 
during  the  transition  from  subsonic  speeds  to  supersonic  speeds.  Also,  ground  combat  vehicles 
can  experience  various  types  of  external  excitations  from  explosives,  gun  fire,  and 
electromagnetic  pulses,  etc.  This  generates  the  need  to  understand  these  types  of  excitations  and 
their  effect  from  a  structural  response  standpoint  for  design  purposes.  Some  of  the  more  common 
transverse  pressure  pulses  are  presented  below. 

The  sonic  boom  over  pressure  which  commonly  materializes  from  exceeding  the  speed  of 
sound  is  expressed  mathematically  as  [5] 

Pt(t)  =  \  (33) 

0  for  f  <  0  and  t  >  rt  ^ 

Where  Pq  denotes  the  peak  reflected  pressure  above  ambient,  t  ^  denotes  the  positive  time 

duration,  from  the  time  of  impact,  and  r  denotes  the  shock  pulse  length  factor.  The  expression  for 

the  sonic  boom  can  be  expressed  functionally  in  contrast  to  a  piecewise  formulation  as 

P,{t)  =  PQ(\-tltp)[H{t)-H{t-rt  p)]  (34) 

Where  H  (t)  denotes  the  Heaviside  Step  Function.  This  implies  that  the  normalized  pressure 
(0  becomes 
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Pmn  (0  =  - (l  -  tjtp  )[h  it)  -  H  {t  -  rt  p)]  (35) 

^mn 

As  special  cases  of  the  sonic  boom,  when  r  =  \ ,  the  waveform  degenerates  into  a  triangular 
pulse.  When  r  =  2,  The  sonic  boom  waveform  becomes  symmetric.  When  r  ^  2  a  non- 
symmetric  N-Pulse  shape  waveform  is  produced.  Finally,  when  r  =  1  and  ->  co  a  step  pulse  is 

realized.  For  the  case  of  the  step  pulse,  the  normalized  pressure  {t)  is  expressed  as 


Pmn  (0  = 


16  Pn 


mn  '  '  2 

mn  K  1 , 


for  t  >  0 


(36) 


As  an  additional  case,  a  rectangular  pulse  can  be  formulated  as 


Pmn  (0  = 


16  Pn 


mn  •'  '  2 

mn  K  1  ^ 


(37) 


Also  under  consideration,  is  a  air  blast  wave  traveling  tangentially/parallel  to  the  plate  which  is 


given  as 


Pj(f)  =  PqC  -  xi) 


(38) 


The  normalized  pressure  P^„  (?)  then  becomes 


Pmn  (0  = 


16  Pn 


mn  K  I  y 


(39) 


Where  c  is  the  wave  speed  within  the  surrounding  medium  which  in  this  is  air;  while  tj  is  a  blast 


decay  parameter.  Finally  a  sine  pulse  is  considered  which  is  expressed  mathematically  as 

i Pq  sm(  ;r  t jt p),  for  0  <  t  <  t p 


Pt(t)  = 


0 


for  t  >  t , 


(40) 
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Functionally  it  is  expressed  as 

Pj{t)  =  Pq  sm(  n:  -  H{t-tp)\ 

Where  the  normalized  pressure  becomes 

16  Pq  /  r  1 

Pmn  (0  =  - Y - tltp)[H{t)-H{t-tp)\ 

^mn 

8.  Results  and  Discussion 


(41) 


(42) 


To  validate  the  present  theory,  validations  are  made  by  simplifying  the  governing  response 
equation  (Eq.  26)  for  the  case  of  an  isotropic  plate  based  on  the  linear  classical  plate  theory 
(CPT)  for  thin  plates.  Both  the  deflection-time  response  and  the  natural  frequency  response  are 
validated.  The  deflection-time  response  for  the  linear  CPT  is  compared  with  Akay[l]  who 
presents  results  for  both  the  linear  and  nonlinear  cases  for  a  constant  uniform  loading.  The 
natural  frequencies  for  the  first  ten  modes  are  compared  against  the  exact  values  provided  by 
Ramu  and  Mohanty  [4]  for  the  linear  CPT. 

8.1  Validation  of  the  deflection-time  response 
As  an  aside,  for  the  case  of  a  uniform  loading  the  transverse  pressure  is  expressed  as 

P  =  Pq  sin(  TTxj Li )  sin(  Tiy  j L2 )  ■  (43) 

In  this  case,  Eq.  (26)  becomes, 

I  w  ('f'l  +  r  w  (t^+K  w  (t)  =  — - — - — 1441 

^  mn  mn  v*  /  '  mn  mn  ^  mn  mn  ^  2  ’ 

71 

While  for  the  case  of  a  isotropic  plate. 
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^11  (x,  y)  =  D  = 


Eh' 


(45) 


12(l-v") 


5£>ii(x,j)  dD^^{x,y)  d^D^^{x,y)  d^D^^{x,y)  d^D^^{x,y) 


dx 


dy 


dx‘ 


dxdy 


=  0 


(46) 


ay- 


d{x  +  /u y,  )  Z.iZ.2  phL^Lj^  . 

=  0  (No  damping  assumed). 


(47a-c) 


Dividing  through  by  gives 


(0 


^mn^mn  C)  P  > 


(48) 


Where  =  \  D [4^  +  Mn)  I P  =  16 Pq  / ph  .  Applying  the  Runge-Kutta  method 

Gives  the  deflection-time  response  which  is  shown  below  in  Fig.  5  as  a  validation  comparing 
Akay  [1]  with  the  present  simplified  isotropic  case.  The  properties  utilized  for  this  validation  are 
given  as, 

Lj  =  L2  =  0.25  m  ,  h  =  0.05  m 

E  =  206  AGPa  ,  p  =  7848  kg  /  ,  v  =  0.25 

Pq  =  0.981  MPfl 

By  observation,  it  can  be  seen  that  very  good  agreement  exists  between  Akay[l],  who  applied  a 
finite  element  technique,  and  the  present  result,  based  on  the  Runge-Kutta  technique.  Any 
differences  are  due  to  the  approximations  made  within  the  two  different  numerical  approaches. 

8.2  Validation  of  the  Natural  frequency 
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To  further  substantiate  the  validation,  comparisons  are  made  with  Ramu  and  Mohanty  [4] 
where  the  first  ten  natural  frequency  modes  for  a  thin  rectangular  plate  based  on  the  linear  CPT 
with  simply  supported  plates  for  various  plate  thicknesses  are  presented.  These  comparisons  are 
listed  below  in  table  2.  The  geometrical  and  material  properties  associated  with  this  data  is  given 
as 

L[  =  0.6m  ,  L2  =  0.4m 
E  =  IQGPa  ,  p  =  2700  kg  /  ,  v  =  0.3 

Very  good  agreement  is  apparent  for  all  ten  modes  for  the  first  two  plate  thicknesses 
(0.00625m  and  0.0125m).  As  the  plate  thickness  becomes  larger,  while  at  the  higher  frequency 
modes,  the  marginal  difference  becomes  larger.  This  is  explained  by  the  fact  that  the  present 
theory  for  thin  plates  is  less  accurate  for  thicker  plates.  Although  the  differences  remain  small, 
the  accurate  prediction  of  the  natural  frequencies  requires  some  of  the  elements  of  the  thicker 
plate  theory  such  as  the  inclusion  of  the  transverse  shear  strain,  normal  strain,  and  the  rotational 
strains. 

8.3  Results  for  the  Present  Theory 

The  present  theoretical  results  are  provided  for  a  variety  of  different  types  of  pressure  pulses 
which  include  a  sonic-boom,  a  triangular  pulse,  a  rectangular  pulse,  a  step  pulse,  a  sine  pulse, 
and  a  traveling  tangential  shock  wave.  The  effect  of  these  pulse  forms  on  the  deflection-time 
response  and  the  effect  of  varying  the  volume  fraction  indexes  is  presented  and  analyzed  to 
evaluate  how  all  of  these  effects  interplay  on  the  structural  response.  The  material  and 
geometrical  properties  adopted  for  the  theoretical  results  are 
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Lj  =  L2  =  0.25  m,  h  =  0.009525  m 

Ti-6A1-4V:  £  =  105  .7  x  10^  Pa,  v=  0.298,  p  =  4429yt^/m^ 

Aluminum  Oxide:  £  =  320  .2  x  10^  Pa,  v  =  0.260,  p  =  3750  kg  /  . 

In  Fig.  6,  The  deflection-time  response  due  to  a  sonic  boom  is  presented  for  three  gradient 
measures  according  to  the  volume  fraction  indexes.  The  case  for  metal  can  be  seen  to  produce 
the  highest  deflections  with  the  lowest  frequency  of  oscillation.  As  the  gradient  distribution 
becomes  more  refined,  the  peak  deflection  values  are  lower  than  its  isotropic  metal  counterpart 
coupled  with  higher  oscillation  frequencies. 

A  traveling  shock  wave  parallel  to  the  plate  and  its  effect  on  the  deflection-time  response  is 
provided  in  Fig.  7.  There  exist  very  high  central  deformations  for  metal  coupled  with  lower 
frequencies  in  contrast  to  ceramic  with  low  deformation  and  higher  frequencies.  Between  the  two 
extremes  lies  the  structural  performance  for  the  multidirectional  case.  Clearly  if  the  desired  effect 
is  to  leverage  the  benefits  of  a  both  material  properties  of  the  constituent  materials,  then 
multidirectional  materials  offer  a  great  advantage  over  isotropic  materials  while  offering  the 
strength  close  to  ceramic.  Another  observation  in  Fig.  7  reveals  that  the  intensity  of  the 
shockwave  decays  much  quicker  for  both  ceramic  as  well  as  for  the  multidirectional  case. 

In  Fig.  8,  the  effect  of  a  degenerated  sonic  boom  into  a  triangular  pulse  in  contrast  to  the 
symmetric  sonic  boom  appears  to  be  less  violent  and  decays  much  faster.  Also  it  should  be  noted 
that  a  comparison  between  multidirectional  grading  and  a  symmetric  unidirectional  grading  can 
be  made.  The  case  [M  =  l,N^  =  N 2  =  0.05  ]  is  very  close  to  a  symmetric  unidirectional 
functional  grading  which  occurs  when  the  volume  fraction  indexes  take  on  the  values 
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[M  =  1,  =  N  2  =  0] .  This  obviously  reveals  that  a  symmetric  unidirectional  functional  grading 

performs  better  from  a  dynamic  response  standpoint  than  its  multidirectional  counterpart. 

In  Figs.  (9)-(l  1),  the  effect  a  step  pulse,  a  sine  pulse,  and  a  rectangular  pulse  are  depicted.  The 
same  kind  of  competing  behavior  is  observed  between  metal,  ceramic,  and  multidirectional 
functional-type  grading.  Fig.  12  provides  a  comparison  between  sonic  boom  loadings  for  various 
magnitudes  of  the  shock  pulse  factor.  When  the  shock  pulse  factor,  r=l,  a  triangular  pulse 
presents  itself,  when  r=2,  a  typical  sonic  boom-type  loading  is  present,  and  when  r=3,  a 
nonsymmetrical  sonic  boom  type  loading  develops.  Comparisons  between  the  three  types  gives 
in  increasing  order  of  shock  pulse  factor  magnitude  the  order  of  most  hostile  response.  The  least 
damaging  is  the  triangular  pulse  with  the  most  violent  response  coming  from  the  unsymmetrical 
sonic  boom. 

Conclusions 

The  dynamic  response  of  multidirectional  simply  supported  plates  has  been  addressed  through 
the  application  of  several  pressure  pulses.  Several  validations  have  been  made  for  the  isotropic 
case.  The  implications  of  the  volume  fraction  index  on  the  material  grading  have  been  presented. 
Also,  various  types  of  pressure  pulses  and  their  effects  on  the  deflection-time  response  of  a  thin 
rectangular  plate  based  on  the  linear  strain  theory  have  been  provided  graphically  and  discussed. 
A  combination  of  semi- analytical  and  numerical  procedures  were  adopted  to  arrive  at  the 
governing  response  such  as  the  Galerkin  method.  The  Gaussian-Quadrature  method  for 
numerical  integration,  and  the  Rung-Kutta  method  to  solve  the  governing  response  equation  for 
the  dynamic  response. 
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It  was  revealed  that  multidirectional  functionally  graded  plates  provide  better  performance  than 
their  isotropic  counterpart  but  decreased  performance  when  compared  with  symmetric 
unidirectional  functionally  graded  plates.  It  is  hoped  that  this  instrumental  base  in  which  to  build 
upon  will  be  realized  in  the  development  of  more  advanced  theories. 
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Figure  Captions 


Fig  1 .  A  pictorial  representation  of  a  simply  supported  multidirectional  functionally  graded  plate 
exposed  to  a  spherical  blast  above  the  plate. 

Fig  2.  The  volume  fraction  as  a  function  of  the  distance  along  the  plate  in  the  x-direction  for 
various  magnitudes  of  the  volume  fraction  indexes  , 

Fig  3.  The  volume  fraction  as  a  function  of  the  distance  along  the  plate  in  the  y-direction  for 
various  magnitudes  of  the  volume  fraction  indexes  N  2  ■ 

Fig  4.  The  volume  fraction  as  a  function  of  the  through-the-thickness  z-coordinate  direction  for 
various  volume  fraction  indexes 

Fig  5.  Validation  of  the  amplitude  of  center  deflection  as  a  function  of  time  for  an  isotropic 
plate. 

Fig  6.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  sonic  boom  for  a 
multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume  fraction  indexes  ( 

=  0,  fp  =  0.005  ,r  =  2,  Pq  =  O.lSMPa  ,  =  0.05  ) 

Fig  7.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  traveling  tangential 
shockwave  for  a  multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume 
fraction  indexes  (rj  =  0.5,  c  =  100  m  /  s,  Pq  =  0.15 MPa  ,  =  0.05  ) 

Fig  8.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  triangular  blast  for  a 
multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume  fraction  indexes  ( 

=  0,tp  =  0.005,  r  =  1,Pq  =  0.15  MPa  ,  A^  =  0.05  ) 

Fig  9.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  step  pulse  for  a 
multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume  fraction  indexes  ( 

=  0,  fp  =  00,  r  =  1,  Pq  =  O.lSMPa  ,  A^  =  0.05  ) 

Fig  10.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  sine  pulse  for  a 
multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume  fraction  indexes  ( 
tp  =  0.005  ,/^0  =  0.15 MPa  ,  A^  =  0.05  ) 
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Fig  1 1 .  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  rectangular  pulse  for  a 
multidirectional  functionally  graded  plate  with  varying  degrees  of  the  volume  fraction  indexes  ( 
tp  =  0.005  ,Po  =  O.lSMPfl  ,  All  =  0-05  ) 

Fig.  12.  Amplitude  of  center  deflection  as  a  function  of  time  due  to  a  sonic  boom  with  varying 
degrees  of  the  shock  pulse  factor  r  for  a  multidirectional  functionally  graded  plate 
(t^  =0,tp  =  0.005  ,  Pq  =  O.lSMPa  ,  An  =  0.05) 
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Tables 


Table  1.  The  effect  of  the  volume  fraction  indexes  on  the  grading-type. 


Grading-Type 

Ai 

Ni 

M 

Volume  Fraction,  Vj 

Isotropic  Metal,  P2 

0(Ai  >  3) 

0  (A2  >  3  ) 

M  >20 

0 

Isotropic  Ceramic,  Pj 

0 

0 

0 

1 

Unidirectional  Symmetric 

0 

0 

1 

Tl  =V^iz) 

Bi-Directional  (x-y  Plane) 

0  <  Ai  <  3 

0  <  A2  <  3 

0 

Tl  =Li(x,v) 

Bi-Directional  {x-z  Plane) 

0  <  Ai  <  3 

0 

0  <  M  <20 

Tl  =V^(x,z) 

Bi-Directional  (y-z  Plane) 

0 

0  <  A2  <  3 

0  <  M  <20 

Tl  =V^(y,z) 

Multi-Directional 

0  <  Aj  <  3 

0  <  A2  <  3 

0  <  M  <20 

Vi  =  Vi(x,  y,z) 
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Table  2.  C  omparison  of  the  natural  frequency  co  (Hz)  for  a  simply  supported  thin  Isotropic 
plate. 


Plate  Thickness 

0.00625m 

Plate  Thickness 

0.0125m 

Plate  Thickness 

0.025m 

Plate  Thickness 

0.05m 

Mode 

Present 

Exact[4] 

Present 

Exact[4] 

Present 

Exact[4] 

Present 

Exact[4] 

m=  1  ,n=  1 

136.4 

136.5 

272.8 

273.1 

545.7 

546.2 

1091.4 

1092.5 

m  =  2,n=  1 

262.4 

262.6 

524.7 

525.2 

1049.4 

1050.4 

2098.9 

2100.9 

m=  1  ,n=2 

419.8 

420.1 

839.5 

840.3 

1679.1 

1680.7 

3358.2 

3361.5 

m  =  3,n=  1 

472.2 

472.7 

944.5 

945.4 

1889.0 

1890.8 

3777.9 

3781.7 

m  =  2,n=2 

545.7 

546.2 

1091.4 

1092.5 

2182.8 

2185.0 

4365.6 

4370.1 

m  =  3,n=2 

755.6 

756.35 

1511.2 

1512.7 

3022.3 

3025.4 

6044.7 

6050.8 

II 

c 

A 

II 

a 

766.1 

766.85 

1532.2 

1533.7 

3064.3 

3067.4 

6128.7 

6134.8 

m=  1  ,n=3 

892.0 

892.9 

1784.0 

1785.8 

3568.1 

3571.6 

7136.1 

7143.3 

m  =  2,n=3 

1017.9 

1018.9 

2035.9 

2037.9 

4071.8 

4075.9 

8143.6 

8151.82 

(N 

II 

C 

II 

a 

1049.4 

1050.4 

2098.9 

2100.9 

4197.7 

4201.9 

8395.4 

8403.9 
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